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Abstract 

A theorem about local in time existence of spacelike foliations with prescribed 
mean curvature in cosmological spacetimes will be proved. The time function of the 
foliation is geometrically defined and fixes the diffeomorphism invariance inherent in 
general foliations of spacetimes. Moreover, in contrast to the situation of the more 
special constant mean curvature foliations, which play an important role in global 
analysis of spacetimes, this theorem overcomes the existence problem arising from 
topological restrictions for surfaces of constant mean curvature. 
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1 Introduction 



Spacelike foliations arise in many contexts in General Relativity. This is related to the 
fact, that the world around us seems three dimensional and we are used to think about 
dynamics as a sequence of processes parametrized by some time function. But in Relativity 
there is no canonical time coordinate, thus any time function only fixes (a part of) the 
diffeomorphism invariance of the theory but this gauge fixing is a priori arbitrary unless 
the time function is tied to some geometrically defined quantity. For that reason, an 
arbitrary time function does not provide any information about the global structure of 
the spacetime under consideration. For example, any thin, spacelike strip in Minkowski 
spacetime can be regarded as global in time, by stretching the time coordinate with an 
appropriate diffeomorphism. This problem is connected with the lack of a canonical metric 
background structure for Einstein's field equations. 

In this work I will define and prove a uniqueness and local in time existence theorem for a 
spatially global, geometrically defined time coordinate in cosmological spacetimes in terms 
of a Prescribed Mean Curvature (PMC) foliation, defined as follows: The leaves {St} are 
determined implicitly by the requirement, that given an initial Cauchy surface Sq := T, 
subject to a certain condition, the mean curvature at each point pt € St is defined by the 
relation H{pt) = H{pq) + t, where pt and po € Sq are connected by the integral curves 
7(t) of the normal vector field of the leaves. Therefore, the elements in this construction 
are purely geometrical in nature, the time function obtained by this procedure depends 
only on the choice of the initial Cauchy surface and we get a geometrical measure of the 
(timelike) size of its Cauchy development as long as this time coordinate exists. I will 
refer to the leaves of such a foliation as surfaces of prescribed mean curvature. Note, that 
there is another definition for surfaces with prescribed mean curvature in the literature 
('H-surfaces'), compare |Q. But unlike the definition given there, where the spacetime 
and a function H prescribing the mean curvature are given quantities, the definition here 
relates to the PMC foliation as a whole, providing an intrinsic prescription without having 
given the spacetime or a function H in advance: Coupling the PMC equations to Einstein's 
evolution equations with matter (such that the Cauchy problem in harmonic coordinates 
is well posed), one gets a system of equations, which constructs a spacetime foliated in 
this way, where the mean curvature 'develops' from leaf to leaf. 

Given this kind of foliation globally, we have a tool at hand for the global analysis of 
spacetimes in an invariant manner, in terms of an asymptotic analysis of a system of 
partial differential equations, throwing away the differential geometric difficulties arising 
from the diffeomorphism invariance. For an introduction to global issues in Relativity 



compare the survey article |An]. Related topics are the development of singularities, the 
cosmic censorhip conjecture and the closed universe recollapse conjecture. Introductions to 
these topics can be found in ||ME[ (global existence and cosmic censorship) and |BT|,|BGT] 



(closed universes and recollapse). 

A related construction are the well-known constant mean curvature (CMC) foliations, see 
the survey article ||^ or | [MT| for reference. Obviously they are a special case of a PMC 
foliation, when the initial Cauchy surface itself has constant mean curvature, and this 
fact sheds some light about the origins of PMC foliations. Indeed, the need for such a 
generalization arises from the fact, that there are strong topological restrictions concerning 
the existence of even a single CMC hypersurface. There exist spacetimes, which do not 



possess any constant mean curvature hypersurface at all, see |Ba]. Conequently, many 



of the results obtained by constructions with CMC foliations presuppose the existence of 
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at least a single CMC hypersurface, which then has to be verified in each concrete case. 
Contrary to this, for the existence of a local in time PMC foliation there is a simple and 
easily verified criterion (compare theorem leading to a broad class of spacetimes, 
possessing at least a local in time PMC foliation. The similarity between PMC and CMC 
foliations then gives rise to the hope, that the global results obtained in the CMC cases 
(under the additional existence assumptions mentioned above) can also be proved using 
PMC foliations, an issue a forthcoming paper will be concerned with. 

The organisation of this work is as follows. The section |2| fixes notation and introduces 
the basic equations and results used throughout this work. Its first part is concerned 
with General Relativity, while the second part is devoted to partial differential equations. 
Section ^ then contains the construction of the local in time PMC foliation. In the last 
section the main result will be stated and discussed. 
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2 Preliminaries 



2.1 Spacetimes and foliations 

A spacetime is a pair {M,g), where M denotes a four dimensional smooth and orientable 
Lorentz manifold with metric g and signature ( — h ++). The metric induces structures 
like the Levi-Civita connection "V and the curvature on M. The sign convention for the 
curvature is fixed by the definition ■'R{X,Y)Z := ^Vx^^yZ - '"Vy^VxZ - ^V^xy\Z, 
where X^Y Z are vector fields. The curvature tensor is then defined as ^K(yV^ X, y, Z) := 
g{W,'^R{X,Y)Z) with Ricci tensor '^Raf3 = scalar curvature ^i? = ''i?^, written 

in abstract index notation of the Ricci calculus. The Einstein tensor reads Ga/d = '^Rai3 ~ 
\'^Rgap and we can write Einstein's field equation as 

(la) Gap = SnTa^ 



or equivalently (by *R = — SvrtrT) 

(lb) X/3 = Svr {Tap - i(tr T) g^p) , 

where Ta/3 denotes the energy momentum tensor of the matter fields. It is a symmetric 
tensor on M with vanishing divergence as a consequence of the Bianchi identities, a re- 
quirement which imposes supplementary conditions on the matter fields coupled to the 
field equation. 

The components of the energy momentum tensor with respect to an observer, represented 
by a unit timelike vector n, have physical meaning: Denote by hap '■= Qap + naUp the 
orthogonal projector on {n}^ in covariant notation, then we can define the energy density, 
momentum density and the stress tensor by 

(2a) p := T^un^v!" 
(2b) jp := -T^.n^'h^p 
(2c) 5„;3 := V/i^/i^ . 



In this work we confine ourselves to cosmological solutions (M, g) of Einstein's field equa- 



tions. Due to |Ba] these are globally hyperbolic and spatially compact spacetimes, where 
the Ricci tensor contracted twice with any timelike vector is non-negative (timelike conver- 
gence condition). This last condition can be reexpressed in terms of the matter variables 
as p + tr S >0 for any observer, known as the strong energy condition. 
An introduction to General Relativity with a treatment of the Cauchy problem for the 



field equations is [Wl], while a recent presentation with a deeper analysis can be found in 



m 



Now let us pay attention to an additional structure. A foliation {St}, t G / C M (/ interval 
containing zero) of (a part of) {M,g) by spacelike hypersurfaces induces on each leaf the 
unit normal vector field n, the metric hap = gap + nanp, which also serves as orthogonal 
projection and the second fundamental form kap := —hah'p'^'S/^ny (the definition of kap 
fixes the sign conventions used in this work). The second fundamental form is a sym- 
metric tensor, intrinsic to the leaves of the foliation, and it can also be written as the 
Lie derivative of the 3-metric h with respect to the normal vector field, kap = —\Cnhap- 
The 3-metric determines further geometrical objects on the leaves, such as the Levi-Civita 
connection V and the curvature tensor R{-). Tensors intrinsic to the leaves of the foliation 
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will carry Latin indices in the abstract index notation. 

The parameter t of the foliation has timelike gradient and thus can be regarded as 
(coordinate-) time. Given local coordinates (x*) on Sq, we can Lie-transport them to 
neighbouring leaves along an arbitrary family of transversal curves, parametrized by x. 
We will express equations containing coordinate components with respect to the adapted 
coordinates {t, x). 

We define the lapse function N and the shift vector i/ _L n on the leaves by the formula 

N = ~g(dt,n) 

(3) dt = Nn + v , thus ^ ' 

u = dt — Nn 

Then we find 1 = dt{dt) = N dt{n). Further, dt is (co-)orthogonal on the leaves and if we 
denote the conormal of the leaves by a we see that dt = —N~^a or a = —Ndt. Thus N~'^ 
measures the length of dt, which can be interpreted as follows: For an observer along n N 
measures the elapsed proper time from St to St+at or equivalently, the elapsed coordinate 
time dt along a journey along n between two leaves, separated by a unit distance (measured 
along n) is A^-\ thus A^-^ measures the elapsed coordinate time along n. Note that these 
relations have only infinitesimal meaning and for the elapsed proper time along n between 
St-^ and St2 we have J*'^ N. 

The coordinates adapted to the foliation together with lapse and shift fix the diffeomor- 
phism invar iance and we get a 3+1-split of the field equations consisting of constraint 
equations intrinsic to the leaves and evolution equations. 

The constraint equations are a consequence of the equations of Gauss and Codazzi, where 
Einstein's equations have been used to eliminate the curvature terms of the 4-geometry. 
The result is 

(4a) R + — = 167r/9 (Hamiltonian constraint) 
(4b) V-'kij — ViH = Sirji (momentum constraint) , 

with = kajjk'^^ and H = \,ik denotes the mean curvature of the leaves. 

The remaining evolution equations for the components of the first and second fundamental 

forms are known as ADM equations. They read in adapted coordinates 

(5a) dthij = -2Nkij + ViUj + VjUi 

(5b) dthj = -ViVjN + N {Rij + Hkij - IklKj - 87r( Sij + \{p-ii S)hij )) 
fk>ij ~\~ kfj \/ j^u ~\~ k'ij- \/ ji/ , 

where the first equation is merely a rewriting of the definition of the second fundamental 
form, k = —^jCnh- Taking the trace of the second equation and eliminating the scalar 
curvature R by the Hamiltonian constraint we obtain the lapse equation 

(6) AN + N(^\kf + 47r{p + tiS)^={dt-iy)H , 

which serves as a constraint of the foliation, induced by Einstein's equations. Note, that 
in cosmological spacetimes, the term in brackets is always non-negative. 
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2.2 Functional analysis and partial differential equations 



At first I state some basic functional analytic definitions. Making use of the multi index 
notation, the Sobolev spaces are defined as 

VF'='P(M") := {/ e LP(M") I a"/ G LP(R")V„|a| < k} , 



and for s € M we define 

H'{W) := {/ G 5'(M") I (1 + l^l'^y/'^Tf G L2(M")} , 

where S' denotes the space of tempered distributions and T is the Fourier transform. 
There are continuous embeddings H^{M.'^) C (M") for s' < s, thus the elements of 
H^{W^), s > are functions. Furthermore, to compare different Sobolev spaces I cite 
[|Ra{ , (Appendix A) where some basic facts from the theory of interpolation spaces are 
stated: 

For So 7^ si in ]R let se := (1 - 0)so + Osi, < 6 < 1. Then 



H 



H 



thus is an intermediate space between i7*o(M") and 

For the convenience of the reader I now list some (well known) inequalities, that are central 
to the analysis of partial differential equations. I adopt the convention of denoting any 
generic constant by C. 

The perhaps most important estimate, we will make use of very often is Gronwall's in- 
equality. 

2.1. Proposition (Gronwall's inequality). 

Let I be an interval, to (z I and a,(3,u G C(/,M+), with 

u{t) < a{t) + / I5{s)u{s)ds 



to 



for all t ^ I. 
Then 



to 



u{t) < a{t) + 
holds for all t £ I. 

The proof of Gronwall's inequality in this particular form can be found in |^ . 

The second important estimate, relating differentiability classes to each other is the col- 
lection of the two Sobolev inequalities, which can be found in standard textbooks about 
partial differential equations or functional analysis, see for example | RR |. Here, ab- 
breviates the vector with components Z?" for all a with s = \a\. 

2.2. Proposition (Sobolev inequalities). 

(SI) l<kp<n 



(S2) kp>n 

In particular we will use a special case of the second inequality, 

2.3. Corollary (Sobolev embedding theorem). 

> § + / 

For the analysis of non-linear equations the Moser estimates are essential 

2.4. Proposition (Moser inequalities). 

(Ml) f,g€ i7*(M") nL°°(R"), s = \a\ 

Wifg^L^ <Cm\L^\\D'g\\L2 + \\D'f\\L49\\L^) 
(M2) / e H'{W) n VFi'~(M"), g £ H'~'^{W) n L~(IR"), s = |a| 

Wifg) - fD^gU. < C {WfWML^ + \\Df\\L^\\D-'^-^g\\L2) 
(M3) / G F^(M") n L~(]R"), F G C°°(M), F(0) = 0, Br := {x £W \ \\x\\ < R} 

\\D'F{f)h.<c{\\f\\L^)\\D'f\\L^ , 
with ciWfU^) = C\\F\\cs(B,,,^^) 



(see for example |Ra| ). 

Let us turn now to the analysis of special classes of partial differential equations. We first 
consider hyperbolic equations: 

2.5. Definition (Symmetric hyperbolic system). 

Let {0} C / C M denote an interval and G C M'^ an open set. Further, for i = 1, . . . ,n 
let A^jA^jB denote smooth functions on / x x G, where — id, A^ and B have 
compact support. A^, A^ are assumed to take values in M.^^^ and B in R^. Then we 
call the system 

, X, u) dtu + A^ {t, X, u) diU + B{t, x,u) = 

for the unknown u : I xG — > R'^ ( quasilinear) symmetric hyperbolic, if the matrices 

A^ are symmetric and A^ is uniformly positive definite. 
The system is called ( inhomogeneous ) linear symmetric hyperbolic, if it has the form 

A°(t, x) dtu + A\t, x) diU + B{t, x)u = f(t, x) . 



The following proposition holds (see e.g. |Ra| for a proof): 
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2.6. Proposition (Existence and Uniqueness of solutions). 



1. Let uq G Co°(M"). Then there exists a unique solution u G C°° ([0,oo[ xM") of the 
linear symmetric hyperbolic system, with u{Q,x) = uq{x). Moreover, u satisfies the 
energy estimate 

\Ht)\\l^<c(\\umls+ sup \\f{t')\\lA 

\ 0<t'<t J 

2. Let s > ^ + 1, G -ff'^(M") with compact support. Then there is aT > which de- 
pends on \\uo\\h' and a unique classical solution u G C"'^([0,T] xM") of the quasilinear 
symmetric hyperbolic system with u{0,x) = uq{x) and 



nG Pi C'' ([0,T],H'-''{W)'^ 

N3k<s 

Moreover the energy estimate 



c {I + \\uit')\\c.) dt' 



1 + Mt)\\h < (1 + 11^(0)11?^.) e 
holds, providing some extension criterion if \\u{t)\\(ji remains bounded in [0, T]. 

In contrast to hyperbolic equations the study of elliptic operators differs substantially from 
that of hyperbolic differential operators, because there is no finite propagation speed, not 
even propagation at all. Therefore we cannot localize and have to analyse elliptic equations 
as a whole instead. 

In this work we need some standard facts about scalar linear elliptic equations on compact 



manifolds. Basics about elliptic operators can be found for example in [RR| and some 
more advanced results on Fredholm theory of elliptic operators on compact manifolds are 
presented in Cantor's article [Q. In particular we find for the Laplacian A := — /i^-' VjVj 
of a smooth compact Riemannian manifold (S, /i): 

2.7. Proposition. 

Let k > be some integer. A : H'''^'^(T,) — > H^^iT,) is a Fredholm operator with 
ind(A) = 0, H'^ = ker(A) Im(A), where the sum is L^ -orthogonal and ker(A) = 
{constant functions}, thus Im(A) = {/ G C°°($]) | / = 0}. 

Moreover, duality and interpolation shows, that A extends to a Fredholm operator 
ff*+^(S) — > H^iTi) of index zero for arbitrary s G M 

The class of Fredholm operators is closed under addition of compact operators and the 
index remains constant. Thus we easily get the 

2.8. Corollary. 

// A G C°°(S) is everywhere non-negative and does not vanish identically, then the 
elliptic operator 

L = A + A : H'+^{^) — ^ H'{J:) 

is an isomorphism for all s G M. 

In particular, this corollary applies to the operator in (0) in cosmological spacetimes. 
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3 Local PMC foliations 



3.1 The PMC equations 

Let (M, g) be a smooth, globally hyperbolic spacetime with compact Cauchy surface S. 
We denote the coordinates on S by and define U to be a Gaussian neighbourhood 
of the form I x S D ] — e, e[ x S, with some e > of S := {t = 0}. Then we have 
g{do,do) = —l,do J- di, i = 1,2,3 and do := dt- D abbreviates the vector (di) and for an 
intrinsic spacelike 3- vector v set := Qahv'^v^ ■ 

Let G J C M denote an interval. A foliation {Sr} of a neighbourhood f/' C f7 of S 
in M by spacelike hypersurfaces is given by a map 99 : J — > C°°(S,]R), r 1 — > ipr, with 
\DipT\ < 1, T I — > (/3r(2;) strictly monotone for all x G S and Sr = {{t,x) \ ^Pt{x) — t = 0}. 
The foliation defines a unit normal n, the first and second fundamental forms h and k 
respectively and the mean curvature H on each leaf, parametrized by r. The coordinates 
on each leaf will be identified with the coordinates on S by the action of with 4>r{x) = 
{(prix),x), so we can regard functions on the leaves as functions on S, parametrized by r. 
With these conventions the foliation fulfills the equation 

(7a) ^ = {l-\D^\^YN , 

where we dropped the subscript r from (/?, now simply considered as a function of (t, x). 
We fix the foliation by imposing a condition for the lapse function A^, which has to obey 

LrN := (A^ + Xr) N = 1 , 
(7b) A, := -/i"^V,V^ , 

Xr ■■= kapk""^ + ^7r{p + trS) , 

which is the usual lapse equation (P) on Sr (indicated by the subscript r in the equations 
above), together with the prescribed mean curvature condition Nn"daH = 1. Thus we 
force the mean curvature to grow along the normal vector field of the leaves as described 
in the introduction. We define a PMC foliation to be a solution of equations ([7a|), (|7b| ) 
and our goal is to prove uniqueness and local in time existence of solutions of that system 
of equations. 

To cast the hyperbolic part of the equations into first order symmetric hyperbolic form, 
we set w := (") := and express all functions in terms of w and N . Then one finds 
n = n{w), h = h{w), k = k{w, Dw) and the equations read 

(8a) drW + A'(w,N)d^w + B{w,N,DN)=0 
(8b) Lriw,Dw)N = 1 , 

with symmetric 4x4 matrices and B takes values in W^. 

Initial values are given by a function wq € -ff*(S,R^), with \vo\ < 1, which represents 
a spacelike hypersurface of regularity in M, with an appropriate choice of s. Since 
we want to have So = S we set wq = throughout this work, but we do not make use 
of this explicit setting in the sequel. Rather we proceed in full generality assuming only 
Wo € F*(S,R^), spacelike. 

In order to get a well-posed initial value problem for the system (^), we have to ensure, 
that Lo(wo, Dwq) is an isomorphism of Sobolev spaces, so that ( pb] ) has a unique solution 
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A^o on '"^0- In view of corollary we require Aq > 0, Aq 7^ on wq which turns out 
to be satisfied for almost all Cauchy surfaces S in spacetimes obeying the strong energy 
condition p + tv S > 0, while A = would force S to be maximal and time symmetric. In 
the latter case a small deformation of T, would suffice. Therefore, from now on we focus 



on cosmological spacetimes in the sense described in section 2.1 



The domain G of definition for the coefficients of the equations can be formally decomposed 
as G = Gw X Gdw X Gat x Gon, where the factors are constraint by the requirements 
\v\ < 1 in and > in R. 

To prove local in time existence and uniqueness for solutions of the system (^) , equivalent 
to (§), we invoke the following procedure: We define iteratively a sequence of uncoupled 
linear equations (see @), which have smooth solutions {w\N^) on non-empty time in- 
tervals . Then we have to show convergence (w^jN^) — > {w,N) in some appropriate 
Sobolev space on a non-empty time interval J. Finally we have to verify, that {w, N) 
represents indeed an unique solution of the system (|8|). 

Let us start with the definition of the sequence mentioned above. Denote by (wq) a 
sequence of smooth functions on S, with \vq\ < 1 and Wq — > wq in H^{T,). The iteration 
will be defined inductively: 



0) Set 



W^{t,x) ■.= Wq{x) y(r,x)(iIxT. 

:= 1 on / X S , 



then one has on := /: \v^\ = \vq\ < 1, which establishes as a well-defined 
family of spacelike hypersurfaces in U. Moreover, A^'' > by construction and 
/ A°(u'°,L»tt;°) > 0, by / Xo{wo,Dwo) > and continuity. 

j) Now we have a non-empty time interval J^~^ for which a family of spacelike hypersur- 
faces w^~^ with 7^ A-' ~^ > is defined, and A^-'"^ > is given on J^~^. Then one 
gets smooth solutions {w^ ,N^) on J^~^ for the (linear) system 

drw^ + A\w^-'^,N^-^) diw^ + B{w^-\N^-\DN^-^) = , 

(9a) 

(0) = 



(9b) Lr{w^-\Dw^-')N^ = 1 ,V,ej.-i , 

since the above conditions on w^~^, X^^^ and A^-^~^ ensure the regularity of the 

equations and the bijectivity of Lt-{w^^^ , Dw^^^). 

Set 

Jj ■= {t e JJ-i I \v^\ < 1, 7^ A^' > 0, > 0} , 

then 2 {0} holds: The compactness of S and |^;q| < 1, 7^ Aq > give l^;-'! < 1, 
7^ A-' > on a non-empty interval by continuity and A^-' > follows from the 
minimum principle for —Lr{w^~^,Dw^~^)N^ = —1. Hence we have constructed the 
desired sequence of solutions {w^ , N^) G G on JK 
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3.2 The energy estimate 



Now we are going to investigate the convergence properties of {w^ , ) . Choose R = 
{Ri,R2), Ri,R2 > appropriate, such that 

Gw := Bii^(lmwo) CC Gw 
Gdw ■■= BR^{lmDwo) CC Gdw 

holds. Moreover, there are domains {1} C Gn CC G^, {0} C Gdn CC Gdn, which 
depend on R, WwqWh"- They will be specified later. 

After this remarks let us turn to the fundamental energy estimate, which will set us in 
the position to prove local existence later on (in subsection 3.3). 

3.1. Lemma. 

Let s > ^ + 2. 

Then there is a T > 0, with J := [—T, T] C , Vj, and a K > 0, such that for all j 
and all r € J: 

\\w^{t)\\hs<K 
\\N^{t)\\hs+^<K 



(b) Im {w^{t),Dw^{t),N^{t),DN^{t)) C x Gd^ x x Gdn • 
Remark. 

For the reader familiar with energy estimates for quasilinear symmetric hyperbolic 
systems I will explain the connection to the standard proofs (see [Ra]). For the 



hyperbolic part of the equations most elements of the standard proof carry over to 
the present situation. The difference is due to the coupling to the elliptic equation 
through the dependence on {N^^^ , DN^^^) in (pa|) and {w^^^ , Dw^~^) in (|9b|). As- 
suming s > I + 2, instead of demanding only s > | + 1 decouples the estimates and 
the proof splits into the standard proof for the hyperbolic part and in a part for the 
elliptic equation. 

Proof (induction over j). 

For j = set T := r° := i diam J°. 
Since w^{t) = Wq we have 

\\w°{t) - wollci = ll^^o - woWci < G\\w^ - woWh^ < Ge 



by assumption and the Sobolev embedding theorem 2.3. Without loss of generality we can 
choose the approximating sequence of the initial values, such that the following conditions 
hold: 

^1 . fRi R2 
^^2"^^Hc^'^, 

WoWhs < 2\\wo\\hs and \\wq^^ - wlWjjs < ^ Vj 
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(we will use the two latter conditions for the local existence proof in Statement (b) 
of the lemma follows then immediately with the remark = 1, DN^ = 0. 
Next, we have ||iV*^(r)||j:^i,+i = = ||l||i2 < C, so, if we set 

K := := max ( ||tt;o||_f/s + e, ||1|[^2 ) 

then (a) follows and we have completed the proof for j = 0. 



Now we want to perform the induction from j — 1 to j. This proceeds in four steps. 
Let f,K be defined ioi {w^ ,N^), . . . , {w^-'^ ,Nj-^). 

Step 1 [Estimate of \\w^ {t)\\h^) 

If a denotes a multiindex of order |a| = s, then differentiation of ( pa] ) yields 

V ' 

I 

^ V ' 

II 

+ {B{w^'\N^-\DN^-^)) . 

^ V ' 

III 

Now we have to estimate the marked terms. Note, that for |r| < T by induction hypothesis 
all arguments lie in a set whose compact closure is contained in G. Furthermore, the 
arguments are bounded in the H^, resp. H^^^ norms, independent of r and j. Then the 
Moser estimates (prop. p.4| ) can be applied and one gets for each |t[ < T 

I <D'^w\l>= I < D'^w^ , A' diD'^w^ >= I < D''w^,{diA')D''w^ > 



< \\DA\\l^\\w^\\]js < C\ 



w 



,i||2. 



\11\\l2 < C(||Z)M||i2||W||Loo + \\DA\\l^\\D''-^Dw^\\l2 
<C7(l + i^)(||w;-''||c;i + ||w;-''||^^.) 

\m\\L2 <cii + K) , 



drh 



hence 



Jt. Jt. 
< C(l + k)[l + \\w^\\ci + Ww^Wh" ) \\w^\\h^ ■ 

Integration yields 

\\w'iT)\\h < \\wlfHs+C{l+K) / {l + \\w^{t)\\ci + \\w^mH^)\\w^t)\\Hsdt 

Jo 
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This estimate remains true, if we replace |r| by T, HtWoHl^s by 2||?i;o|l//s, (t) ||ci 
C\\w^ {t)\\H'' (Sobolev embedding theorem, corollary |2.3| ) and terms of the form (1 + a)^ 
with C(l + a^), so we get 

\\w^iT)\\ls < 2\\wofHs +C{l + k) r {1 + \Wmls ) dt . 

Jo 

Gronwall (prop. |2.1[ ) then gives the final estimate 

\\w^{T)\\l.<{l + 2\\wo\\h)e''^'^^^^ , V,-V|,|<^ . 

Set now 

:=max(KO,(l + 2||«;o|||,.)e^) 
1 



Tw ■= min ( T 



^0 



then we have < K^, < and for all j and |r| < \\w^{t)\\h'' is bounded by K^, 
hence we have statement (a) of the lemma with respect to w. 

Step 2 {Boundedness of \\w^ {t)\\qi) 

The equation ( pa| ) gives an Sobolev estimate for drW^: 

< \\A'diW^\\Hs-i + \\B\\hs-i < C{1 + K^f , 

since < C{l+Kyj)Kw and < C{l+Kw). Further, with multiindices 

/3,7, = s — 1, I7I = s — 2 the differentiated equation yields 

\\D'<drDw^\\L2 = WdrD^W^U^ < ||I||i2 + ||II||i2 + ||III||i2 , 

with a replaced by P in each of the terms I, II, III defined above. Since we have s > | + 2, 
we can estimate II and III as before with s replaced by s — 1, and for I we have 



||I||i2 = [ < A^diD^w^ ,A'diD^w^ >< C\\A\\loo\\w^\Is < CK 
and together 

\\drDw^\\Hs-2<CKl + C{l + K^){l + K^)<C{l + K^f 

It follows 

\\w\t) - Wo\\l°^ < \\wl - WoWl^ + I \\drW\t)\\Loc dt 



2 

w 1 



^0 

1 







< ^uAn{Ri,R2) + C{l + K^f\T\ 

\\Dw^{t) - DwqWl-- < \\Dw^Q - DwqWl-o + / \\drDw^{t)\\L^ dt 

Jo 

<^mm{Ri,R2) + C{l + K^f\T\ . 



If we set 



. / - imin(i?i,ii2) 
T,„ := mm T,„, — 
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then we have < and for all |r| < ||^^-' (''") — '"^ollci is bounded by mm{Ri, R2), 
which means, that statement (b) of the lemma is fulfillled with respect to w. 

Step 3 {Estimate of \\N^ {t)\\hs+i) 

We are looking for a bound of \\L-r{w^~^ , Dw^~^)~^\\ in the H^~^ — H^^^ operator norm. 
Lt-{w^~^, Dw^^^) is a sum of terms of the form ai~^D'^, with < 2 and at~^ = 
ar{w^~^ , Dw^~^) where Ur is a smooth function of its arguments and the parameter r. 

We proceed in several steps. 

At first we consider Lt-{w^~^ ,Dw^~^) as an operator i7^(E) — > H^{T) = L'^{T>) only. We 
already know ||Lo(it''', -Du;'^)"^ || < C and show successively 

1) \\Lo{w^,Dw^)-^\\ < C Vj 

2) \\Lr{w^ T Dw^)^^\\ < C Vj V|t-|<7^j^, with some T/v > chosen appropriately. 

The proof of these statements relies on the following standard fact: 
Note. 

If T : — > i?0(S) is invertible and S : H^{T.) — > H^{T.) satisfies 



IT -511 < iir 



-iii-i 



Then S = T - {T - S) = T{1 - T-^{T - S)) is invertible with norm \\S-'^\\ < 
\\T-'\\ Er=o^' = T^,\\T~'l where q = WT-^T - S)\\ < 1. 

Proof of 1): 

Lq{'w, Dw)~^ exists by assumption. Now define 5jLo := Lo{w,Dw) — Lq{w^ ,Dw^). 
We want to show |j5jLo|| < e for an < e < \\Lq{w, Dw)~'^\\~^ . The note then 
completes the claim. 

Since Wq — > wq in we can choose the sequence such that for the coefficients 

of (5jLo ll'^ifloll/fs-i < e holds (where the Moser estimates (see prop. 2A) have 
been applied on some expression, coming from the mean value theorem, compare 
the note in subsection |3.3| . In particular we have for / G ff^(S) the estimate 
PjLo/llj^o < C||5jao||L°° ||/||_H"2 < and we are done. 

Proof of 2): 



We have shown in 1), that \\Lq(w^ , Dw^) ^\\ < C, hence C ^ < \\Lq{w^ D 
We set 5t-L{w^ ,Dw^) := Lq{w^ ,Dw^) — Lt-{w^ , Dw^) and must show, that the esti- 
mate \\6rL{w^,Dw^)\\ < C-^ holds. 

For / G H^iT.) and |k| < 2 we have \\6ra^ D'^ fW^jo < C||(5T-a-' ||l°° ll/ll_f/2i and it is 
enough to find a bound for H^T-a-' • For this we mention, that 6ra^ = dra^ ■ 
Thus 

||5^a-'||ioo < llmaxl^T-a-' I T||ioo = max(max|(9T-a-' [) T 

r XT 

= max ( max I 5t- a-' I) T = maxH^T-a-' T < C{1 + Kyj)T 

T X T 

Finally choose some 0<Tn <f with C(l + K^)Tn < . 

Until now we have proven ||A^-^(t)||j:/2 < C\Lr(w'-'^^ ^Dw^^^)^^\ < C, Vj, V|^|<j'j^, where 
Lr{w^^^ , Dw^^^) acts as an operator H'^{Ti) — > H^{Ti). What remains to do is to get 
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control in higher norms. For this note, that the coefficients ai ^ of Lt{w^ ^, Dw^ ^) have 
regularity |[ai~^||j:/s-i < C{\+K) and therefore < C{\+K). Now let us consider 

Lr{w^~^,Dw^-^){DN^) = D[Lr{w^-\Dw^-^)N^] - {DLriw^-\ Dw^~^))N^ 

= -{DLr{w^-^,Dw^-^))N^ 

and \\{Dai-^)D^N^{T)\\L2 < \\Dai~^\\L--\\N^ {t)\\h2 < C{1 + K). Therefore, the right- 
hand side of the equation above lies in -L^(S) and is bounded there independently of j 
and r, Lr{w^-'^ , Dw^''^) is an isomorphism H^{Y.) — > ^^(S), thus DN^{t) G H'^{T.), 
bounded by 

\\DN\t)\\h2 < \\Lr{w^-^,Dw^-^)-^\\\\{DLr{w^-^,Dw^-^))N^\\L2 <C{l + k) 

thus ||-/V-'('7")||h3 < C(l + K), and we have control over N^{t) with one differentiability 
order increased. 

Now we iterate this procedure to get estimates on higher norms. For multiindices /3, 7, 
with I7I = — 1 we find 

= D{Lriw^'\ Dw^-^)iD^' N^)} - {DLr{w^-\Dw^-^))D''N^ 

K+A=7 

where successively in each step the previous results has been inserted. Since we have a 
bound at most for ||ai~^ we get the highest order estimate on = s — 1. Thus 

we have D'-^N^{t) G H^{^), with \\N^t)\\hs+i < C{1 + K). Setting Kn := C(l + K) 
completes the proof of (a) of the lemma with respect to A'^. 

Step 4 {Boundedness 0/ ||A^-' (r)||ci) 

Multiplication with the Sobolev constant (compare cor. 2.3) for H^^^ ^ yields 



\\N^{t)\\ci<CKn , 

= :Rn 

SO we can set 

Gn := Br^{{1}) , Gdn ■■= BR^m) . 

These domains are independent of r and j and well defined, N^{t), DN^{t) are contained 
in Gtv, Gdn for all j,|T| < T, giving statement (b) of the lemma with respect to N and 
setting 

T := min(r^,rAr) < f 
K := max(i^tt,, K^) 

finishes the proof of the lemma. ■ 

Inspection of the proof shows, that we have shown more. In particular we have the 

3.2. Corollary. 

3 
2 

W{t)\\c^.<Rn , 
hence lTa.{D'^N^) C Gn, for \a\ < 3 and all j and \t\ < T. 



With s> ^ + 2 and \\N^ {t)\\hs+i < it follows 
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and 

3.3. Corollary. 

For all j and \t\ < T the operator Lr{w^ , Dw^) : H'^{T,) — > L'^(J^) is invertible and 
the estimate 

\\Lr{w\Dw^)-^\\ < C 

holds. 

Moreover, the constants T and K depend only on R and HwoIIj?"- 
3.3 Local existence 

The aim of the present subsection is to establish the local in time existence theorem, by 
proving convergence of the sequence {w^ , A^-' ) towards a solution of the system (^) . Since 
S is compact, we will get spatially global results. 

With the notation 

6jU := u^~^^ — 

for any given sequence {u^) we have to estimate the differences 6jW and 6jN in Banach 
spaces of the form 

(c'{J,X),\\-\y^j)^x) , := J^supp^txllx , J C J := [-f,f] , 

where T denotes the time constant of the energy estimate and X an appropriate Sobolev 
space. Cauchy sequences with respect to this topology are uniformly convergent in r. 

As a consequence of the mean value theorem we 

Note. 

For U CR'',V CR^ let F :U xV — ^ a C^-function. 

Then there is a continuous function M on (U x V)"^, with M ee eM„ e W^i^+i)^ 
M„ G M™^^'^, M„ G M.'^''^, such that for all ui,U2 G U, vi,V2 G V 

F{U2,V2) - F{ui,Vi) = M ~ ""M = M.,{U2 - ui) + Mu{v2 - vi) 

\V2 -VlJ 

holds. 

This kind of replacement will be used frequently in the sequel. 

After these preparations the equations for the differences will be derived: Fix r G J, then 
it follows from = drW^ + A\w^-^ ,N^-^) diW^ + B{w^-^ ^N^'^ .DN^-^) = drW^"^^ + 

= dr5jw + A\w^,N^) didjw + A\w^,w^-\n\N^-^) diW^ ( ^^^In ) 
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hence the equation fulfilled by the difference djW is 

= drSjW + A\w^,W) di6jw 

(10a) / Sj^iw \ 

+ B(w^,Dw^,w^-\N^,DN^,N^-\DN^-^){ 5,-i7V 

\5,_iDnJ 

and i? is a smooth function of its arguments. 

For N one gets in the same manner from 1 = Lt-{w^~^ , Dw^~^)N^ = Lr{w^ , Dw'-')N^^^\ 

= Lr{w^ ,Dw^)N^^^ - Lr{w^-\Dw^-^)N^ 
= Lr{w^ , Dw^)N^+^-Lr{w\ Dw^)N^ +Lr{w^ ,Dw^)N^ -Lr{w^-^,Dw^-^)N^ 

= Lr{w\Dw^)6jN - friw^,Dw\w^-\Dw^-\N^,DN^,D^N^) (s%\'dw) > 
hence the equation for 5jN is 

(10b) Lr{w\DuP)5jN = fr{uP , Dw^ ,up-\ Dw^-\ , DN^ , D^N^) ^^p^^^ ) 

and Jt is a smooth function of its arguments and its parameter r. 
Inspection of the last equation yields the important 

3.4. Lemma. 

Let s> 1 + 2. 

Then there exists a constant C > 0, such that 

\\5,N{T)\\H2<C\\5,.rw{T)\\H^ , V.V^gj 

holds. 

The proof of this lemma is an easy consequence of the fact, that the right-hand side of 
(|10b| ) is an element of -L^(S) and Lt-{w^ , Dw^) : if^(5]) — > L'^i^) is an isomorphism, with 
bounded inverse by corollary 3.3. 

3.5. Remark. 

With this lemma at hand, we achieve again a decoupling of our estimates. More- 
over, Cauchy sequences {6jw) in C^{J,H^{T,)) imply Cauchy sequences (SjN) in 
{J , H'^ (T,)) . Since ||5jA^(r)||j|^s+i and ||5jt(;(r)||/i-s are uniformly bounded, we 
get the inequality ||5jA''(r)||Hr < C\\6j-iw{T)\\fjr-i uniformly in r G J for all 
l<r<s + lby interpolation. 



Now let us turn to the estimate of the sequence {Sjw) in {J , (T,)) . Applying D 
on (|10a| ) we infer from the energy estimate that all coefficients of the differentiated 
equation are bounded in all relevant norms (independent of j and t £ J). So we get the 
following estimate (in the same way as we have done so far in the energy estimate for 
equation (|9a])). 

— \\D8jw\\\2 = -r < DSjW,DSjW >= 2 / < DSjW,drDSjW > 
dr dr J^, 

< C ( ||D(5jt(;||2,2 + ||Z)(5j_i'u;||/^2 + ||Z)(5j_iA^||j:/i ) ||Z)(5jt(;||£,2 , 
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hence 



-^¥jw\\m < C ( ¥iw\\m + ¥j-M\m + ¥i-iN\\H-^ 

< C ( WdjwWjjr + \\6j-iwfH^ + \\6j.2w\\jji ] 



where lemma p.4| has been used in the last step. Integration yields for all r E J (this 
means |r| < T), using the abbreviation 5jWo '■= Wq'^^ — uPq: 

Jo 

< ( ¥jWo\\Hi+C {\\6j.iw{t)\\jj, + \\6j-2w{t)\\jj,)dt \ e^^ , 



where the Gronwall estimate (prop. |2.1| ) has been applied in the second step. Since all 
operations were independent of r, we can improve the estimate to 



COiJ), 



Choose a r > 0, such that CTe^"^ < ^7 < 1, and set J := [— T, T] C J. Then we get 

ll'^i^lIcO{J),L2 ^ ll^j^ollHi + 9 ( ll'^i-l^llcO{J),Hi + ¥i-M\l'i[j),H^ 
summation over j = 2, . . . , m and rearranging of the terms gives 

m 



J'^\\C0{J),H 

3=2 ^ ... 



^1100(7), HI. 



Clearly, the right-hand side of this inequality is bounded (by choice of the approximat- 
ing sequence of initial values, we have Pjti'oll^i ^ see the j = 0-step in the proof 
of the energy estimate), independent of m, so we have found a majorant for the left- 
hand side. Thus the sequence {w^) is a Cauchy sequence, hence converges in the Banach 
space C^{J, H^iTi)). Moreover, standard arguments for the interpolation of Sobolev norms 
yield convergence in C^{J,H'' {T.)), for all s', with 1 < s' < s, since 5jW is bounded in 
CO(J,^r"(S)): 

w := lim in {J , H'' {T.)) , 

j^oo 

and the Sobolev embedding theorem asserts, that w is a classical function of regularity 
C0(J,C2(S)) c C70(J X S). 

It remains to show, that is a solution of (|8a|) . For this purpose we investigate the 
sequence {dr^jw). Equation ( |10ei| ) yields 



< C ( \\6jw\\cQ(j)^H^ + ||5j-l-w||cO(J),//i + ¥j~lN\\cO(j)^H^ ) 

< C ( ¥jw\\co(_j)^H^ + ||5j_i-u;||co(j),/fi + ||5j-2w||co(j)^Hi ) 
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with C again independent of j and r by the energy estimate. This shows, that {drW^) 
is a Cauchy sequence in C'^(J, ff^(S)) and since drSjW is bounded in C^{J, H^~^(T,)) 
interpolation defines 

w:=limdrW^ in C\j,H''-\^)) , 

and the Sobolev embedding theorem guaranties, that is a classical function, lying in 
CO(J,CnS)) C C70(J X S). 

The convergence of the sequences (w^) and {drW^) is by construnction uniform in r, hence 
we have w = drW for w G C°(J, n C^{J, H'''^{J:)) C C^{J x S). Moreover, a^-u; 

obeys 

d^w = w= lim a^t;;^' = lim (-A\w^-'^,N^-^)diW^ - B{w^-^ ,N^-^ .DN^'H , 
with smooth and therefore C'^( J, i7*'~^(S))-continuous A^,B. 

We conclude, that w is a solution of equation (^a]), if the sequence (A^-^) converges in 
CO(J,i/^'(S)). 



Let us now come to the inspection of the sequence {5jN). Lemma 3^ and remark 3.5 
show Pj-^||(70(j) //s'+i < C!\\5j-iw\\(~,o(^j-^ Hs' and it easily follows 

Af:=limiV^' in C°{J, H''-^\^)) . 

The smoothness of the coefficients in Lt-{w^~^ , Dw^~^) , t ^ J ensures, that for all / G 
j^s +1(^5])^ Lt-{w^~^, Dw^~^) f converges in the norm to Lr{w, Dw)f, which means 

convergence Lr{w^~^ , Dw^'^) — Lt-{w, Dw) in the H^''^^ —H^'~^ operator norm. From 
this we can infer Lt-{u!^~^ , Dw^~^)N^ — > Lt-{w, Dw)N uniform with respect to r in 
since 

\\Lr{w,Dw)N - Lr{w^~^,Dw^-'^)N^\jj,>^-, 

< \\Lr{w,Dw) - Lr{w^-\Dw^-^)\\ ||iV(T)|lj:^,/+i 
+ \\Lr{w^-\Dw^-')\\ \\N{t)-N^{t)\\j,,,+, 

and both summands consist of one bounded and one converging to zero term. The only 
possible limit for Lr{'w, Dw)N is 1, since Lr{'w^~^ , Dw^~^)N^ = 1 for all j. Hence 
solves equation (pb|). 

Furthermore, the Sobolev embedding theorem gives regularity C^(S) for each N{t), and 
since w{t) € C2(S), Dw{t) € C^(E), Lr{w,Dw) acts as an operator C3(S) — > C^C^), 
establishing N{t) as a classical solution of ( |8b| ) for each t £ J. 
Altogether we have proven the following 

3.6. Theorem. 



Using the conventions described in subsection p. J| given the system (|8D with spacelike 
initial values wq G H^(T,), s > |+2 and 7^ Aq > 0, there exist a non-empty interval 
J = [—T,T] and classical solutions {w,N) of @, satisfying 

w G c°{j,H''{j:))ncHJ,H''-\j:)) n c\j x s) 

for each s' G M vuith 1 < s' < s. 
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3.4 Improving the regularity 



Now we want to get rid of the primed s, as well as to investigate the effect of time 
derivatives on w and N, in order to make the results more satisfying. 



Standard arguments on symmetric hyperbolic equations (see |Ra|) show, that w{t) pos- 
sesses indeed regularity H'^, since only the energy estimate and the convergence of the 
sequence {w^) are involved in the arguments. For N{t) we simply replace s' + 1 by s and 
we get 

3.7. Note. 

so we have a well-defined initial value mapping wq i — > {w{t), N{t)), mapping iJ*(E) 
into itself for all t & J. 

We want to know something about the effect of differentiation with respect to time. For 
we have seen, that applying dr on decreases spatial differentiability by one. For A^-' 
we see first, that in H^{T,) we get (compare corollary |3.3| and step 3 in the proof of the 
energy estimate N{t) = limj N^(t) uniformly in r, and since N^{t) G H^^^{T,) this 
also holds in i?** "''^(S) by interpolation. The same conclusion for the sequence {drN^ (t)) 
holds for s' -|- 1 replaced by s', since the coefficients of the operator drLT-{w^~^ , Dw^~^) 
have regularity decreased by one. Thus {drN^{T)) converges uniformly with respect to r 
in (S) and this is enough to show, that drN{T) = limj drN^(T) in (S): 

3.8. Note. 

N G C\j,H''+\^))nC\j,H''{^)) 

Let us combine the previous results. Iteration of the calculations show, that the effect of 
time differentiating w or is a decrease of spatial differentiability order by one and we 
get the 

3.9. Corollary. 

iw,N)e n C'(J,i7^-'(S)) , 

1<S 

in particular 

w;oGC°°(S) =^ {w,N) eC^iJ x^) . 



3.5 Uniqueness 

Finally we consider the question of uniqueness. Let two classical solutions (wajNa), a = 
1, 2 of the system (^) with initial value wq G //*(S) be given, with regularity 

Wa G c\j, H%i:)) n c\j, i/^-i(s)) 
NaeC\j,H''+H^)) 
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for s > I + 2, < s' < s. 

Define fj, := wi — W2, v := Ni — N2, ttien fj, obeys the difference equation 



drfl + A'{wi,Ni)difl + B{Wa,DWa,Na,DNa)[^J =0 . 

This is an inhomogeneous linear symmetric hyperbohc equation of the form drf^ + A'^difj,+ 
Bfj, = f, whose coefficients depend after insertion of the solutions {wa,Na) only on (r, x) 
and f{T,x) = f{v{T,x),Dv{T,x)). The energy estimate for equation of this type yields 
(see prop. ^ ) 



MrW^, < c [wmrHs' + sup 

V 0<t<T 

In the same way v fulfills an analogous difference equation, and can be estimated like 
< C\\L{wi,Dwi) - L{w2,Dw2)\\ < ( < C\\fi\\fjs' 

\\Du\\j^,, <C\\(^Di,^y\jjs'-2 <CMjjs' , 

for each r. This simplifies the energy estimate for to 

Mr)\\l.> <c(mO)\\1,+ sup \Ht)\\l.' 

V 0<t<T 

Since we have ^(0) = by assumption, we find /u(t) = for all t ^ J. Furthermore 

||z^(r)||^a'+i < C\\i-i{t)\\^s' yields 1^ = 0, hence 

Wi = W2 

Ni = N2 , 
as desired. 
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4 Conclusion and outlook 



Putting together theorem |3.6| , corollary 3^ and the uniqueness result we get the final 
4.1. Theorem. 

Let (M, g) he a smooth, globally hyperbolic spacetime, obeying the strong energy 
condition, with compact Cauchy surface S and 

A = +47r(p + tr5) 

does not vanish identically on E. 

Then there exists aT > and a unique smooth PMC foliation {Sr}, t £ [—T,T] in 
{M,g), with T, = Sq. 



Note, that the setting here is quite general, no symmetry assumptions have to be made 
and essentially the strong energy condition turned out to be sufficient for the local in time 
existence of a unique PMC foliation up to the choice of an initial Cauchy surface. 

What remains to do is to globalize the result. Here are two problems involved: How large 
is the interval of values taken by the time coordinate and does the global foliation then 
cover the whole spacetime? To answer these questions in general there seem to be no 
techniques available up to now. One strategy to obtain global results is, to study first 
spacetimes with some spatial symmetry, taking advantage of the simplifications of the 
equations. Then the hope is, that the techniques developed in these cases give insight into 
the nature of more general classes of spacetimes, by successively lowering the degree of 
symmetry. This line of attack will be the content of a forthcoming paper. 
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